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[57] ABSTRACT 

A method of and system for rendering a halftone image 
of a gray scale image by utilizing a pixd-by-pixcl com- 
parison of the gray scale image against a blue noise mask 
is disclosed in which the gray scale image is scanned on 
a pixel-by-pixel basis and compared on a pixel-by-pixel 
basis to an array of corresponding data points contained 
in a blue noise mask stored in a PROM or computer 
memory in order to produce the desired halftoned im- 
age. Both digital and optically implemented halftone 
methods are disclosed. Application specific modifica- 
tions of the blue noise mask as well as its use for produc- 
ing halftoned color images are also disclosed. 

11 Claims, 10 Drawing Sheets 
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METHOD AND APPARATUS FOR HALFTONE 
RENDERING OF A GRAY SCALE IMAGE USING A 
BLUE NOISE MASK 

5 

This application is a continuation-in-part of U.S. pa- 
tent appUcation Sen No. 07/622,056, filed Dec. 4, 1990, 
now U.S. Pat. No. 5,111,310 and entitled METHOD 
AND APPARATUS FOR HALFTONE RENDER- 
ING OF A GRAY SCALE IMAGE USING A BLUE 10 
NOISE MASK. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to the halfton- 
ing of images. More particularly, the present invention 15 
relates to a method of and system for rendering a half- 
tone by utilizing a pixel-by-pixel comparison of the gray, 
scale image against a blue noise mask. 

Many printing devices are not capable of reproducing 
gray scale images because they are bi-level. As a result, 20 
the binary representation of a gray scale image is a 
necessity in a wide range of applications such as laser 
printers, facsimile machines, lithography (newspaper 
printing), liquid crystal displays and plasma panels. 
Gray scale images are typically converted to binary 25 
images using halftone techniques. Halftoning renders 
the illusion of various shades of gray by using only two 
levels^ black and white, and can be implemented either 
digitally (facsimile machines, laser printers) or optically 
(newspaper printing). 30 

Halftoning algorithms are classified into point and 
neighborhood algorithms according to the number of 
points from the input gray scale image required to cal- 
culate one output point in the output binary image. In 
the case of digital halftoning, pomts correspond to pix- 35 
els. In point algorithms, the halftoning is accomplished 
by a simple pointwise comparison of the gray scale 
image against a nonimage, usually aperiodic (but not 
always) array or mask. For every point in the input 
image, depending on which point value (the gray scale 40 
image or the mask) is larger, either a 1 or 0 is placed 
respectively at the corresponding location in the binary 
output image. Halftoning using neighborhood algo- 
rithms is not done by simple pointwise comparison, but 
usually requires fdtering operations that involve a num- 45 
ber of points from the input gray scale image in order to 
calculate one point in the output image. 

At present, given the existing halftoning algorithms, 
the choice for a specific halftoning algorithm depends 
on the target device and always requires a trade-off 50 
between image quality and speed. Neighborhood half- 
toning algorithms result in a good image quality (al- 
though the image is not completely artifact free), but 
they are slow and cannot be optically implemented. 
That leaves point algorithms as the only choice for 55 
optical applications such as newspaper printing. Point 
algorithms are fast and are suitable for all target devices, 
but the output usually suffers from artifacts such as 
periodic artifacts and false contours. 

The halftoning system disclosed herein utilizes a 60 
point algorithm, and combines the output image quality 
of neighborhood algorithms with the speed and wide 
application range of point algorithms. A point algo- 
rithm is utilized and the halftoning is achieved by a 
pixel wise comparison against a nonimage array, called 65 
the "blue noise" mask. 

The digital halftoning of images with multiple levels, 
such as gray scale levels, is known in the art Two major 
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techniques are currently in use. They are the ordered 
dither and the error diffusion methods. See Digital Half- 
toning by R. Ulichney, MIT Press, Cambridge, Mass. 
(1987). See also R. W. Hoyd and L. Steinberg. "Adapt- 
ive algorithm for spatial gray scale", SID International 
Symposium Digest of Technical Papers, pps. 36-37. 
The Floyd and Steinberg paper is directed to the digital 
halftoning of a gray scale. 

The major ordered dither techniques are the clus- 
tered-dot dither and dispersed-dot dither techniques. A 
white noise random dither technique is seldom utilized 
because it produces the poorest quality image and, of 
the other two dither techniques, clustered-dot is by far 
the most used. Both of those techniques are based upon 
a threshold screen pattern that is of a fixed size. For 
example, 6x6 threshold screens may be compared with 
the digital input values. If the input digital value is 
greater than the screen pattern number, a 1 is produced 
and, if it is less, a 0 value is assigned. The number of 
levels that can be represented using either technique 
depends on the size of the screen. For example, a 6x6 
screen can produce 36 unique levels. 

More levels can be achieved with larger patterns, 
however, a reduction in the effective resolution occurs 
because the ability to transition among levels is at a 
coarser pitch. At the pixel rate of about 300 to 500 per 
inch, which is the average pixel rate of copiers and laser 
printers, the pattern artifacts are visible for screen pat- 
terns larger than 4X4, and, smce 16 levels do not pro- 
vide an adequate precision for typical continuous-tone 
imagery, a suboptimal resolution is usually obtained. 

One solution to such a problemi is disclosed by Ulich- 
ney in a paper **Dithering with Blue Noise" pubHshed in 
the Proceedings of the IEEE, Vol. 76, No. 1, January 
1988. In that article, a method of spatial dithering is 
described which renders the illusion of continuous-tone 
pictures on displays that are capable of only producing 
binary picture elements. The method produces a blue 
noise pattern high frequency white noise from a filter to 
provide desirable properties for halftoning. More spe- 
cifically, Ulichney uses perturbed-weight error diffu- 
sion methods which when digitally implemented run at 
a much slower speed (approximately 100 times slower) 
than is attainable with the present invention. 

Error diffusion techniques, such as that disclosed in 
the Ulichney IEEE article, are fundamentally different 
from ordered dither techniques in that there is no fixed 
screen pattern. Rather, a recursive algorithm is used 
that attempts to correct errors made by representing the 
continuous signal by binary values. 

The error diffiision meiod described by Ulichney, 
and others, such as Floyd and Steinberg, also has the 
disadvantage that it requires scanning, convolution- 
style calculations and, although it can be implemented 
for use with copiers, facsimile machines, etc., requires 
local calculations. It cannot, however, be optically im- 
plemented. In addition, all error diffusion techniques, 
including those described by Ulichney and Royd and 
Steinberg, show scanning and start-up artifacts, which 
are not present in the instant invention. Also, while 
Ulichney describes a method that produces blue noise, 
the blue noise patterns produced by the present inven- 
tion are more isotropic than those produced by Ulich- 
ney or other error diffusion methods. Utilizing ordered 
di&er methods produces notably periodic patterns that 
are even much more obtrusive than those produced by 
error diffusion methods. 
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In some prior art systems, all dot profiles correspond- 
ing to different gray levels were derived independently, 
as if each grade level was its own special case. Thus, for 
example, in U.S. Pat. No. 4,920,501, to Sullivan et al., 
many individual dot profiles, corresponding to the de- 5 
sired number of gray levels, must be stored. In the pres- 
ent invention, on the other hand, dot profUes are built 
"on top of the profiles from lower gray levels, such 
that a single valued 2-dimensional function, that is, the 
cimiulative array or blue noise mask, can be con- 10 
structed. When that single valued function is thre- 
sholded at any level, the resulting binary pattern is 
exactly the blue noise dot profile design for that particu- 
lar gray level, p(i j,g), where p can be one or zero corre- 
sponding to black or white, i and j are coordinates of IS 
pixels, and g represents a gray level 0<g< 1. 

Another drawback to prior art methods is that the dot 
profile for a given gray level was designed to have blue 
noise properties by indirect methods, such as using an 
error diffusion filter with perturbed weights (Ulichney) 20 
or by a "simulated annealing" algorithm, as in U.S. Pat 
No. 4,920,501. The method disclosed herein is advanta- 
geous wiih respect to the prior art in that the desired 
blue noise power spectrum is produced through the use 
of a filter on the dot profile and is implemented directly 23 
in the transform domain. Such filtering results in a 
nearly ideal blue noise pattern with implicit long-scale 
periodicity because of the circular convolution implicit 
in the use of discrete Fourier transforms. However, the 
filtered pattern is no longer binary. Thus, a minimiza- 30 
tion of error approach is utilized in which the largest 
differences between the ideal, filtered, blue noise pat- 
tern and the unfiltered dot profile are identified. The 
magnitude and location of those differences indicate the 
pixels in which ones and zeros could be changed to 35 
produce a more ideal blue noise dot profile. 

Display devices, including printing devices as well as 
media, have their own unique input-output characteris- 
tics. In some uses, such as medit^ ultrasound imaging, 
the user has traditionally been provided with some 40 
control as to the final gray scale mapping. For example, 
the user may be able to select between high and low 
contrast images. The display and film characteristics, 
however, must be accounted for in each rendition. 

In the area of halftone rendering, traditional halftone 45 
screens using small (for example, 8x8 pixel kernels) 
provide only limited degrees of freedom to alter the 
input-output characteristics and usually a linear cumula- 
tive distribution function (CDI^ has been reported. See 
Digital Halftoning by R. Ulichney, MIT Press, Cam- 50 
bridge, Mass. (1987). See also, R. Bayer, **An Optimum 
Method for 2 Level Rendition of Continuous Tone 
Pictures", IEEE International Conf. Comm., 1973, and 
G. C. Reid, Postscript Language Program Design (green 
book), Addison-Wesley Publishing Co., New York, 55 
N.Y. (1988), page 137. By a linear CDF, it is meant that 
\0% of the halftone kernel pixel contents will be less 
than 10% of the maximum value and that 50% of the 
pixels will contain values less than 50% of the mayimitTp 
values, and so forth. 60 

In the case of the blue noise mask method disclosed 
herein, a large unstructured pattern of, for example, 
256X256 pixel kernels, provides sufficient degrees of 
freedom to modify the cumulative distribution function 
so as to provide both linear and non-linear mappings of 65 
input and output characteristics. That makes it possible 
to construct specialized blue noise masks in which a 
particular printer output and media characteristics can 
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be compensated for by a modified blue noise mask gen- 
erated as disclosed in this application. 

The present inventive method herdn may also be 
applied to color halftonings by independently threshold- 
ing each of the component colors against the disclosed 
blue noise mask and then overprinting. Such method 
produces a pleasing pattern without any blurring of the 
image. Such method is a great improvement over the 
known prior art, which is discussed below. 

In U.S. Pat No. 5,010,398, there is disclosed a 
method for color corrections by dry dye etching using 
a photographically produced oiask which may be used 
in the production of printing plates for printing repro- 
ductions of colored originals and in which a contact 
print is overexposed to a photographic mask. The pho- 
tographic mask is constituted so as to isolate a selected 
area in addition to being exposed normally for obtaining 
an exact copy of an original halftone separation. The 
mask is electronically generated by scanning each sepa- 
ration, digitizing each signal and then storing the digital 
values in a digital storage device. 

U.S. Pat. No. 4,974,067 relates to a multi-step digital 
color image reproducing method and apparatus which 
separates an original image into a plurality of color 
components to produce image data associated with each 
respective one of the color components. The image data 
are individually processed to provide record color com- 
ponent density data which data are used to record a 
halftone representation pattern of that color compo- 
nent. 

An apparatus and methods for digital halftoning is 
disclosed in U.S. Pat. No. 4,924,301 for producing half- 
tone screens or color separations from continuous tone 
intensity signals that are supplied by an optical scanner. 
Using a digital signal processor, the continuous tone 
intensity values are processed to estabhsh memory maps 
which, in conjunction with a digital data output device 
such as a laser printer, produces the desired halftone 
screen. The digital signal processor utilizes a dither 
matrix in order to produce halftone screens having a 
screen angle that does not substantially differ from the 
screen angles of the yellow, cyan and magenta color 
separations in conventional four color halftone printing. 
A dither array is also utilized to produce the halftone 
screens having a screen angle that substantiaUy differs 
from the screen angle used in the black halftone color 
separation in conventional four color halftone printing. 

U.S. Pat No. 4,342,046 relates to a contact screen for 
making color separation halftone blocks for use in a 
picture reproducing machine in which a plurality of 
halftone screens having different screen angles are ar- 
ranged on a base film in the corresponding positions of 
color separation reproduction pictures to be repro- 
duced on the base film and transparent blank spaces are 
formed between two adjacent halftone screens. 

A method and apparatus for making monochrome 
facsimiles of color images on color displays is disclosed 
in U.S. Pat No. 4,308,533 for making 35 MM color 
shdes from a color image created on a color cathode 
tube terminal. U.S. Pat. No. 3,085,878 relates to the 
preparation of traditional halftone screens for color 
separation. 

U.S. Pat No. 4,657,831 relates to the production of 
electrophotographic color proofs of hal^ne dot pat- 
tern images which closely simulate the dot gain of prints 
made with lithographic plates and liquid inks. 

A process for the production of photographic masks 
is disclosed in U.S. Pat. No. 4,997,733 in which such 
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masks are used for the tonal correction by dry dot etch- It is another object of this invention to provide a 

ing in which the selection of a particular halftone color system for the halftone rendering of a gray scale image 

separation image or overlaying registering combination which tailors the image thus produced to compensate 

of halftone color separation images is determined on the for particular output printer and media characteristics 
basis of optical density differences in at least one such 3 such that the undesired display of media characteristics 

halftone color separation. Such differences include dif- can be substantially eliminated, 

ferences in contrast, between each area to be isolated as It is yet another object of this invention to provide a 

a substantially transparent area, and at least one particu- system for the halftone rendering of continuous tone 

lar background area surrounding each area to be iso- ^o\ot images such that pleasing images are produced 
lated. 10 with little spreading out of the color and with no blur- 

U.S, Pat No. 4»477,833 is directed to a method of ring of the image, 

color conversion with improved interpolation in which Briefly described, these and other objects of the in- 

an apparatus for converting a color image from one vention are accomplished by generating a blue noise 

colored space to another colored space mcludes a mem- which, when thresholded at any gray level g, 

ory which stores a finite number of output signals which produces a blue noise bmary pattern appropriate for 

represent colors in the output space and which is ad- S^ay level. After the blue noise mask has been 

dressed by signals representing a color in the input generated, it is stored in a PROM. The image to be 

space. The interpolation process is utilized in order to halftoned is then read by a scanner on a pixel-by-pixel 

derive an output color value for an input color located compared to the corresponding stored pixel in 

between colors stored in the memory. ^^^^ ^ produce the resulting binary 

A method and apparatus for producing halftone "^fi® converted to a 

printing forms with rotated screens based upon ran- ^P^y ^^^^ ^ resultant halftoned image, 

domly selected screen threshold values is disclosed in ^ generated and 

U.S. Pat No. 4,700,235. The screens have arbitrary '^^^^^^^ ^V^y modified to tailor it for a particular 

screenanglesandscreen width. Screendotsareexposed output printer and inecha cWtem^ 

on a recordmg media by means of a recordmg element Pf^f ^^n is provided m a blue noisemask for undesired 

whose exposure beam is switched on and off by a con- ^^Z^J h it ^l^^^^f ^« ™f *^ 

trol signal modified by altermg the first order statistics or the 

NonVof the foregoing references have the advan- ^X^rnw!"^^^^^ 

tagesoftheuseofthebluenoisemaskmethoddisclosed S^^^^t and media dependent effect. Tlie 

h^ in^oduc^g pleasing isotropic, non-clumpy ^of XT^^ tcT^^T^;^^^^^ 

moire resistant patterns with only some spreadmg out of ^^^t, ru^^^ui^^ +1,^ iLi«^ ^""'^ . ^ r^"^ "*"t ^^r^ 

^^1^^ ;^ir „^*t. ui — ; At. • m2sk. Punchmg the blue noise mask mvolves settmg the 

the color or ink but with no blurnng of the image. ^^^^^ .TLi„.-^ ^ ^ - • • i r 

£, • extreme low values to a certain mitnm^if p value, for 

SUMMARY AND OBJECTS OF THE 35 example, 0 and setting the extreme high values to a 

INVENTION certain maximum value, for example 255. The values 

In view of the foregoing, it should be apparent that ^^^J"' 
there still exists a need in the art for a method of and The meihod of generating and utilizing the blue noise 
apparatus for the halftone rendermg of gray scale mi- 40 mask discussed above can also be ^pU^ to color half- 
ages m which a digital data processor is ^rtihzed m a toning, by independently thresholding each one of the 
simple and precise manner to accomplish the halftone component colors against the blue noise mask and then 
r^dermg to provide a bmary scale image which is overprinting the halftoned component color images, 
characterized by the pixelwise comparison of the image jhe blue noise mask can also be shifted by one pixel 
being analj^^ against a blue noise mask. TJere also 45 before it is used on each of the different color planes. In 
exists a need for the modification of such halftone ren- that manner, the color energy is spread out over a larger 
doing process m order to counter some undesirable space For example, the blue noise mask can be shifted 
prmter and media dependent effects such that the half- one pixel up or to the side, when the red image and blue 
tone rendermg of gray scales is modified to produce image are being halftoned, respectively. Other varia- 
mput and output characteristics unique to a particular 50 tions and modifications for using the blue noise mask for 
type of device or media and in which the images pro- color halftonmg are discussed in the specification. Such 
duced are superior to those produced without such principles may be used for either RGB halftoning or 
modification process. CMYK color printing. 

It should likewise be apparent that there still exists a in an optical implementation, the gray scale image is 
need m the art for a method for color halftoning by 55 photographed through the generated blue noise mask 

independently thresholding each of the component and the resultant superposition is printed onto a high 

colors against a blue noise mask m order to produce a contrast fihn. An additive photographic process may 

pleasing, isotropic, non-clumpy, moire resistant pattern also be utilized in which the blue noise mask is added to 

with only some spreading out of the color or ink but the gray scale image at the film plane, for example, by a 
with no blurring of the image. 60 double exposure process. The photographic blue noise 

More particularly, it is an object of this mvention to mask can be obtained from a calculated blue noise array 

provide a system for the halftone rendering of a gray using a fihn prmter interfaced to the PROM or com- 

scale image which has a simple and reliable mechanism putcr in which the blue noise mask array is stored, 
for producing the desired image. 

Still more particularly, it is an object of this invention 65 DESCRIPTION OF THE DRAWINGS 

to provide a system for the halftone rendering of a gray FIG. 1 is a drawing showing the power spectrum of 

scale image which can be implemented cither digitally a blue noise pattern formed in accordance with the 

or optically. present invention; 
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FIG. 2 is a diagram of a flow chart for the design of 
the blue noise mask of the present invention; 

FIG, 3 is a diagram of a flow chart for the digital 
implementation of halftoning using a blue noise mask in 
accordance with the present invention; 

FIG. 4 is a schematic block diagram of a hardware 
system for digitally implementing halftoning using the 
blue noise mask in accordance with the present inven- 
tion; 

FIG. 5 is a drawing of a multiplicative photographic 
process utilized for optically implementing halftoning 
using a blue noise mask in accordance with the present 
invention; 

FIG. 6 is a drawing of an additive photographic pro- 
cess which may be utilized in the optical implementa- 
tion of halftomng using a blue noise mask in connection 
with the process shown in FIG. 5; 

FIG. 7 is a diagram of a flow chart showing the modi- 
fication of a blue noise mask to produce a punched, 
linearized version of that blue noise mask; 

FIG. 8 is a diagram of a flow chart for the modifica- 
tion of a blue noise mask using the concave down sigma 
curve modification to produce a high resolution version 
of the blue noise masl^ 

FIG. 9 is a diagram of a flow chart for the modifica- 
tion of a blue noise mask using the concave up sigma 
curve modification to produce a low resolution version 
of the blue noise mask; 

FIG. 10 is a drawing showing the relationship of the 
number of pixels to the value of pixels for a linear blue 
noise mask; 

FIG. 11 is a drawing showing the relationship be- 
tween the number of pixels and the value of pixels for a 
non-linear^ high contrast blue noise mask produced after 35 
applying the CDSC direct mapping processing with 
punch; and 

FIG. 12 is a diagram of a flow chart showing the 
application of a blue noise mask to color halftomng. 
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Prior to referring to the drawings, the following de- 
scription of the theoretical underpinnings of the present 
invention is provided. 

As described above, the present invention is a half- 
tone rendering system which accomplishes its function 
by a pixel-by-pixel comparison of a gray scale image 
against a "blue noise" mask. As referred to herein, the 
term **h\vie noise*' is a pattern with negligible low fre- 
quency components which, possesses certain visually 
pleasing properties, as described by R. Ulichney in his 
book, Digitai Halftoning. 

In the present invention, depending upon which pixel 
is larger, either the gray scale image or the blue noise 55 
mask, a 1 or a 0 is placed in the binary (black or white) 
image file which is the halftone rendered version of the 
gray scale image. Using the notation that the gray scale 
unage is M X N pixels in size and B-bits of gray per pixel, 
the blue noise mask can be a smaller array JxK in size 60 
where J is less than or equal to M and K is less than or 
equal to N with only B— 1 bits per pixel. 

The blue noise mask described herem is constructed 
to have unique first and second order properties. When 
thresholded at any level, for example at A% of the 65 
maximum level, exactly A out of every 100 pixels will 
be greater than the threshold value. In addition, the 
spatial distribution of the pixels above the threshold is 



arranged in such a manner as to form a blue noise pat- 
tern which has been shown to be visually pleasing. 

The disclosed blue noise mask, therefore, has the 
characteristic that the first order statistics are uniformly 
distributed over gray levels. That is, when the blue 
noise mask is thresholded at a gray level g, exactly 
gX 100% of all values are below the thre^old. For 
g— 0.5, exactly 50% of the blue noise mask pixels are 
above, and 50% below the threshold value. The blue 
noise mask disclosed herein also has the characteristic 
that when thresholded at any level g, the resulting bit 
pattern has a power spectrum consistent with and ap- 
proximating the ideal blue noise pattern for that thresh- 
old. In addition, since the blue noise image is con- 
structed with explicit 'Svraparound" properties, a small 
blue noise pattern of J X K pixels can be used to halftone 
render a larger MxN pixel's image, because the pixel- 
by-pixel comparison can proceed modulo J and modulo 
K in the respective directions, with no apparent discon- 
tinuities or obvious periodicities. However, the value of 
(JXK) should not be smaller than X/2, where X is the 
number of levels of the original gray scale image. 

It is also desirable to describe the digital haUtoniiig 
system of the present invention for the analog case in 
which discrete space is replaced by continuous space. 
Using such notation, x and y represent continuous 
space, while i and j represent discrete space. Thus, the 
gray scale image is denoted by f(x,y), the blue noise 
mask is denoted by m(x,y) and the output (halftoned) 
binary image is denoted by h(x,y). 

Thus, for a B-bit image array f(ij), the blue noise 
mask array m(i j) is a B-bit array such that, when thre- 
sholded against f(i,j), up to 2^ levels of varying distribu- 
tion of black and white dots can be represented on a 
rectangular grid. Note that the dimensions of the blue 
noise mask can be smaller than those of the gray scale 
image and that the halftoning of the gray scale image is 
achieved by a periodic repetition of m(ij) over the 
entire image plane. For example, for a 256x256 8-bit 
class of images, a 128 X 128 8-bit blue noise mask array 
can be used. 

The binary pattern that results after thresholding the 
blue noise mask at a constant level g is called the dot 
profile for that level. The dot profiles are arrays that 
have the same dimensions as the mask array, and consist 
of ones and zeros. The ratio of ones to zeros is different 
for every dot profile and depends on the gray level that 
particular dot profile represents. In the notation used 
herein, the higher the gray level, the more ones and less 
zeros that will be contamed in the dot profile. p(i j,g) is 
used to denote the value of the dot profile at pixel loca- 
tion (ij) and for the gray level g. g=0 is used to repre- 
sent black and g»l is used to represent white. Thus, 
O^g^l. Also, by denoting as fij the value of the dis- 
crete space function fOJ) at pixd location (iJ), a NxN 
binary image h(x,y) can be written as follows in terms of 
the dot profiles: 



(1) 



nR)/R 



where R is the spacing between the addressable points 
on the display device, and rect(x)=l is |x|<| and 
rect(x)=0 otherwise. Therefore, for any gray scale 
image, the corresponding binary inaage h(x,y) can be 
constructed as follows in terms of the dot profiles: For 
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every pixel in the gray scale image array fCiJ) that is at flow chart showing the steps for designing a blue noise 

the (i j) location and has a value fy— g, the correspond- mask for generating level g+ Ag from level g. 

ing pixel in the binary image array h(ij) has a vahie that Step 1, Take the 2-dimensional Fourier transform of the 

is given by the value of the g-level dot proffle at the (ij) dot profile ppj,i] and obtain the dot profile P[u,v,i], 

location. 5 where u and v are the transformed coordinates, and P 

The dot profiles for every level are designed and represents the Fourier Transform, 
combined in such a way as to build a single valued Step 2. Apply a blue noise filter D(u^v,i) to the spec- 
function, the blue noise mask. The blue noise mask is trum P[u,v,J] and in that way obtain new spec- 
constructed such that when thresholded at any level, trum P'[u,v,J]=P[u,v,i) x D(u,v,i). The blue noise 
the resulting dot profile is a locally aperiodic and iso- ^0 filter is designed to produce in the dot profile speo- 
tropic binary pattern with small low-frequency compo- ^^^^ P'[u,v, J] an average cross section along a radi- 
nents, which in the halftoning literature, is known as a symmetric line shown in FIG. 1. The principal 
blue noise pattern. Those dot profiles are not indepen- frequency is given by 1 V2 R. 
dent of each other, but the dot profile for level gi -h Ag ^- '^^^ Inverse Fourier transform of P'[u,v,i] 
is constructed from the dot profile for level gi by re- obtain p'[ij,i], which is no longer binary but has 
placing some selected zeros with ones. For example, for ^^^^ ^^^^ properties, 
a NxN B-bit mask array and maximum pixel value difference e[iJ,i]-p'[ij,J]-p[ij,i]. 
given by 2^, Ag is given by Ag= and the number of difference is referred to as the error array, 
zeros that wiO change to ones, in order to go from level ^^^P ^* Classify all pixels into two classes according to 
gitolevelgi-|-AgisN2/2^. ^ the value of p[ij,i] for each pixel; all the zeros belong 

As the dot profile is changed from its pattern at gi to ^ ^ second. Then, 
gi -t- Ag, another array called the cumulative array is ™^ ^^^^ ^ ^® ^^^^ ^^^^ ^ classes accord- 
incremented in such a way as to keep track of the c< l^^^ ^P*'**^ P^®^ 
changes in dot profiles from gray level to gray level. ^^^P >w=€^or die magnitude of the high- 
That cumulative array (not a binary array but a B-bit ^ T acceptable error. That Imut is usually set equal to 
array) becomes the blue noise mask because, when thre- f magnitude error. For the zeros. X,=€ 
sholded at any level g, the resulting binary pattern rc ^t^Z ^^^^^ "^-^u ^"f^ 
produces the dot profile for that leveL ^^^T J^. 

Referring now to the figures wherein like reference Jf^^^V^^^^ f "^1 ^^f 

numerals are used throughout, there is shown in FIG. 1 ^t contam a one and have an em^^ 

o Ai^^^ «f fi,^ «^«,«. i» ^uuwu m rivj. x defined negative hnut to zeros. The number of zeros 

a dut^ of the power spec^mi of a blue noise pattern ^^ ^^ ^ ^ ^ ^ 

^ITTIjI^^^^^^^ averageispreserved.TlieinitiaIizationprocessisthen 

ponents m the frequency domam corresponds to the 35 complete 

absence of disturbing artifacts in the spatial domain. ^he abov^e procedure is then repeated until no pixels 

Radi^ymmefry m the frequ^icy domam corresponds w an error higher than some ^determined error, 

to ^tropy m ttie spaUal domam. Isotropy, apcnodiaty Note that the magnitude of the average error becomes 

and the lack of low-frequency artifacts are all desurable lower for both zeros and ones every time the procedure 

properties m halftonmg because they lead to visually 40 is repeated. 

plying patter^. In order to finish the initialization procedure, refer to 

As shown m FIG. 1, the cutoff frequency fg, which is another NxN array, which is denoted as c[ij,i] and 

termed the Prmcipal Frequency, depends as follows on referred to as the cumulative array, and give a value of 

the gray level g: 2-^-1 to every pixel whose corresponding pixel in the 

45 dot profile has a value of zero, and give a value of 

(2) 2^-1—1 otherwise. In that way, when the cumulative 
array, which eventually will become the blue noise 
mask, is thresholded at a 50% gray level, the resulting 
dot profile is equal to p[ij,ij. 

SO Alter having generated in the above fashion the dot 

where R. as before, is the distance between addressable ^^^Z^^l^V'^^^' *^f, ^tt^ ^I^Z^ ^ 

J' t > « . •^ww^ »w«*wi>wx*^ constructed, where Ac is usually taken as i-^. the aufiuiti- 

S^nl^l can be se«i from the above equa- j + Ag gray level is generated from the dot profile for 

tion f^ac^bevesitsnmrnn v^^^ 55 the i level by converting zeros to ones. TTie 

that level ^e populations of black and white dots are selection of the pixels th^ contain a zero and wiU be 

equal ^ thus very high frequency components appear repkced by a one is done foDowing a procedure shnilar 

m tiie J°^ge. , . , , to the one described previously for the design of the i 

For a NXN B-bit image with 2^ as the maxunum dot profile in FIG 2 

pixel value, ^e blue noise mask is constructed as fol- 60 In general, the dot'profile for the g+ Ag level can be 

lows: First, the dot profile p[io. J] that corresponds to generated from the dot profile for the g level, as shown 

the 50% gray level is created. That dot profile is gener- in FIG. 2. Up to Step 4, the procedure for the creation 

atcd from a white noise pattern after filtering it with a of the g+ Ag dot profile is exactly the same as the pro- 

highpassed circularly symmetric filter and results in a cedure for the creation of the initial dot profile for the 

binary pattern having visually annoying low frequency 65 J level. It is important to note that in Step 2, the princi- 

con^onents. In order to give blue noise properties to pal frequency of the blue noise filter is updated for 

the p[i j, J] dot profile, the following iteration procedure every level according to equation (2). After Step 4, the 

is utilized, similar to that shown in FIG. 2, which is a purpose is to go up one gray level and thus only zeros 
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are changed to ones. Using the error array, the pixels than N— 1. That indicates the end of a row or column, 
that contain a zero are classified in Step 5 and rank If a negative determination is made at step 312, then j is 
ordered, and then W'/2^ selected zeros are changed to set equal to j-f- 1 at step 314 and the program then re- 
ones in Step 6: turns to step 304 to scan the next pixel. 

^ r . . , . , . . ^ \f2Ji affirmative determination is made at step 312, 

PCv.5]=on^v^s/rt'IU^+^l=i;^>i 0) -^^^^^ 3^^^ ^ 

Every time a zero is changed to one. the statistics of ^^J^^?/ '^^^^^ instant method is applied to the^ftrst 

its neighborhood change and therefore the information ^=0) of the next hne. A determmation is then 

contained in the error array for its neighboring pixels in ^^e at step 318 of whether i is greater than N- 1. If an 

may not be valid any more. For that reason, only a few affirmative determmaUon is made at step 318, that mdi- 

zeros are replaced with ones and then the error array is ^ ^^^^ reached, and 

recalculated or an additional criteria is checked, such program then ends at 320. 

as, ndghborhood mean and nmlengths. Finally, the a negative determination is made at step 318, that 

cumulative array is updated in Step 7 by adding one indicates that the end of the image may not have been 

only to those pixels that still correspond to a zero in the reached and that additional pixels remain. Thus, the 

dot profile p[iJ>g+Ag]: next line is scanned. The value of i is then set equal to 

i+1 at step 322, the value j is set equal to zero at step 

c[y.^+A«J=«[Ud+;»lw-ff+A^J-- 8>i (4) 316 and then the next pixel is scanned at step 304. 

,Q FIG. 4 shows an example of the hardware which may 

where the bar indicates a logical •'not" operation chang- -'^ used for the digital implementation of halftoning using a 

ittg zeros to ones and vice verea. bluenoisemaskasdiscussedintheinstantappUcation.lt 

In that fashion, when the blue noise mask is thre- should be understood that the hardware implementation 

sholded at a constant evel g+ Ag, theresultmg binary ^^^^ ^ ^1 ^ ^ 

pattemis the dot profUe p[iO,g+ Ag]. That procedure is operational amplifier in place of the comparator 402 in 

repeated mitd Ae dot proffl^ for dl the gray levels ^ ^ ^ significant advantege of digital 

from i-hAg up to 1 are created. The levels from J — Ag u ■ • ui • i *i. * •* • i. r ^ 

to 0 are created in the same way with the only differ- ^^r^^""! ' blue noise mask that it is much faster 

ence that the ones are changed to zeros and the cumula- ^ ^^^."^ blue noise producmg techmques, 

tive array is updated as follows: ^^^^ halftonmg is done by a simple pixelwise 

30 comparison. The digital halftoning using a blue noise 

As] =<iiJ»g\-plU'g~^l' g< I (5) mask process of the present invention can be speeded up 

even further by performing the comparison step in par- 

When the process has been implemented for all gray aUel, since all of the thresholds are preassigned. 

levels g, the cumulative array contains the desired blue As shown in FIG. 4, a scanner 400 is used to scan an 

noise dot profile for all levels, and is therefore the de- 35 image and convert the pixels on that image from an 

sired blue noise mask. array of f(x,y) to the gray scale image array f(ij). The 

Once the blue noise mask has been generated, as de- output from die scanner 400 is fed to a first input of the 

scribed in connection with FIGS. 1 and 2, it can be used comparator 402. 

in a halftoning process. Since halftoning using a blue As previously described, the blue noise mask array 

noise mask is a point algorithm, it can therefore be im- 40 mfi j) is stored in the PROM 406 which may be located 

plemrated either digitally or optically. a computer 405. The output from the PROM 406 is 

A flow chart depicting the digital nnplementation of fed to the second input of the comparator 402. The 

halftoning using a blue noise mask according to the ^^^^^ ^he comparator is the binary image array h(i j) 

present mventonis shown m FIG. 3. In digital apphca- ^^^^ ^f^^ 404 which converts that 

tions such as facsimile machines and laser prmteis, the 45 j^e final image array h(x,y). 

instant method requires much less memory and/or com- ^ previously discussed, the present halftoning using 

putation than do other blue-nois^p^^^^ T^^'iv a blue noise m^k invention can also be implemented k 

such ^^OTor diffusion tau^^^^ an optical or photographic mamier An e^ple of an 

van et al m U.S. Pat No. 4,920,501, issued Apr. 24, T^^i v puuiutjiapmw ujoaacr^ /ui cjuuopjc ui an 

. ' opUcal apphcation of the present halftoning system is 

Thenecessarymemoryneededtostorethebluenoise Photographic process used in new2)aper printing, 

mask array is stored on a PROM, as shown in FIG. 4. ^^f^^ t^'^T' ^ multiphcative or additive. 

Then, the halftoning of a N X N gray scale image array ^ '"'i^P'f?^? photographic process, the gray 

f(i j) against the blue noise mask array m(io) is imple- "^Se K^fV) 500 is photographed through the blue 

mented as follows: The i and j variables are first set to 0 55 ^^^^ ^ ^ transmittance m(x,y) and 

at steps 300 and 302, respectively. The next pixel fy is ^® resultant superposition h(x,y)=f(x,y) x m(x,y) is 

then scanned at step 304. A determination is made at printed onto high contrast fflm 504, such as high gamma 

step 306 to determine if the value of that pixel fy is procedure is shown in FIG. 5. It should be 

greater than the value of the corresponding element m understood that a point ^ in the array f(x,y) corresponds 

the blue mask array m/y. 60 to a dot h;, in the array h(x,y), whose size and shape 

If it is determined at step 306 that the value of the depends on the gray level that f), represents, 

gray scale image array pixel is greater than the value of FIG, 6 shows the additive photographic process in 

the blue noise mask array pixel, then the value of the which the blue noise mask is added to the gray scale 

resulting array h/^ is set equal to 0 at step 310. image at the film plane by a double exposure. The gray 

If an af^rmative determination is made at step 306, 65 scale image array and blue noise mask array are added 

then the value of the resulting array element hij is set by adder 600 and then fed to the high gamma film 504 

equal to 1 at step 308. After steps 308 and 310, a deter- which produces the halftoned output. The adder 600 is 

mination is then made at step 312 of whether j is greater in simplest form a repeated exposure of the film 504, 
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where the image and the blue noise mask are exposed Ught and dark regions. That will result in the produc- 

separately onto the fibn 504, which is then developed. tion of a higher contrast image. 

In general, a ^yscklephotogr^hic blue noise mask Themapping function operates on pixels of the blue 

m(x,y) ^ be obtamed from a calcukted array m(ij) noise mask, b(ij), that arc equal to a specific value, g. 

^#niS^.^/^ ""f ^ and sets that specific value iual to a new value, g': 
a PROM or a computer. The conversion from discrete 

points to a continuous image is then given by equation for an 6(W)-«^'(w)=^ (8) 
(1). Print film or transparency film is exposed by the 

computer coafroUed fihn printer so as to produce a where f(g)=g' is a single valued, non-linear mapping 
photoffra^^^ 10 fmiction chosen to alter the image rendering, and b^? 

T^eblueno^emaskcanak^ is the output, modified blue ndse mask. 

^^.^^^f "^"^f"'' ^ "^1?" FIGS. 7-9 show diagrams of the flow charts for pro- 

The digital halftomng process for a bmary display using a.^pJIZ *7 c • 

a blue noise mask ^G. 2) can also be pressed J ducmg three different modifications of ^e blue noise 

follows; <^ »^ cApioacu jj^y^ ^jjg p^^jj Yg^^Qjj^ jjjgjj contrast 

15 concave down sigma curve (CDSC) and the low con- 

A(g)=toKm(^+Av)} (6) concave up sigma curve (CUSQ. All versions 

include the **punch". Punching the blue noise mask 

wherein int denotes integer truncation and the gray means that the extreme low values are set to a certain 

levels of m(ij) and f(ij) vary between 0 and 1. In gen- minimum value, such as Q, and that the extreme high 

eral, for a K-bit display, the output image array h^J) values are set to a certain mgifiTniiTn value, such as 255, 

can be written as follows: and the values between the maximum and mir^i Hi HT n 

values are then re-linearized. 

h^L A - 1 i^ffnK n„n ft _l ^^ U) ^® CDSC modification of the blue noise mask is 

- IFTT"^^ - «Mw) + Aw)> accomphshed by setting F(g)=g3 which produces a 

'Tiigh contrast*' mapping. The CUSC modification sets 

The 2^— 1 threshold values are equally spaced between f(g)=to g-^, thus producing a flatter or low contrast 

0 and 1. A non-uniform quantizer is also possible. curve. 

It is also possible to modify the blue noise mask in After the modification of the blue noise magir is ac- 

order to minimize undesirable printer and media depen- complished, halftoning of the desired image using that 

dent effects. That may be accomplished by modifying mask is done in the manner previously described. That 

the first order statistics or the cumulative distribution is, the halftoning is performed by comparing the values 

fimction (CDF) of the blue noise mask (BNM). Such of the image with the values of the modified mask. If the 

modificatoon is useful in an environment such as medical value of the image is larger, the value of the new image 

ultrasound mia^g m which the user, using the present is rendered as black. Otherwise, the value of the new 

mvention, may be able to se ect b^een high and low 3^ ^ ^ ^^^^ ^ ^^^^^^ 

contrast images and the display and fihn characteristics a more compact form (reduced from 8-bits) because 

of such medical ultrasound miagmg eqmpment can be there is only one bit per pixel in a binary or halftone 

accountedformeachrendition.Inthec^ofbluenoise j^^e. In the modified masks, each S k SiTSue 

n^ a large unstructured pattern of, for example ^obe mask is mapped to a new value based i^on Se 

256x256 pixel kernels provides sufiBaent degrees of ^ „a^«ui^ v "'^w jaiuc ^^w" luc 

freedom with which to modify the CDF so as to pr^ ^ '^T"'' "^^^ Punch^DSC and 

vide both linear and non-linear mappings of the iiput Halftomng is Aen p^ormed based on those 

and output. TMt makes it possible to^M spe*^- ''Z l'lr^' f T"?^ * values, us«s can 

ized blue noise masks forVparticular output printer. ^lodify the miage to ehmmate artifacts and other defects 

Medium characteristics can also be minimized in such 45 ^« ^P^^ Pnnter and/or media being uti- 

modified blue noise masks. Therefore, as will be de- , 

scribed hereafter, the present invention provides for the halftone screens created using the modified 

altering of the cumulative distribution function of the characteristics of the blue noise 

blue noise mask that would otherwise be utilized so as to non-linear input-output characteristics, 

produce unique and more desirable input-output char- 50 While the gray scale characteristics of the blue noise 

acteristics. While three such examples of such modified mask are modified, the isotropic, unstructured, visually 

blue noise masks are discussed herein, those of ordinary pleasing pattern of black and white pixels produced by 

skill in the art will readily recognize other ways in such blue noise masks are preserved, 

which to modify the blue noise mask in order to achieve ^ * diagram of the flow chart for modifying a 

similar results. 55 blue noise mask to produce a linearized version of the 

The first order statistics of the blue noise mask can be blue noise mask. At step 700, the value b(i j)=g of 

obtained directiy from the mask itself as previously the blue noise mask to be modified is read. At 702, the 

discussed. The method for modifying the mask involves maximum and mmitri^iTrf values specified by the user are 

taking each individual value of the mask and mapping it obtained. The values of the blue noise mask obtained at 

into a new value (known as direct value mapping), 60 step 700 are then re-lineaiized at step 704 by setting 

while avoiding certain extreme points of the blue noise r =255/(max-min) and then calculating the values of the 

mask. In that manner, the samc^ non-clumpy image can modified blue noise mask b'(ij)=((b(ij))— min)*r. In 

be produced. The mapping can be performed in such a this example, we assume an 8-bit blue noise mask is 

manner that some undesirable printer output character- being used, hence, max and min wiU lie between 0 and 

istics can be eliminat ed. For instance, when the output 65 255, respectively. 

device produces low contrast knages, the mapping can A determination is then made at step 706 of whether 

be chosen to modify the mask and the CDF of the mask each value of the modified blue noise mask b'(i j) is 

in order to enable more pixels to be depoated in tiie greater tiian the max value. If an affirmative determina- 
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tion is made at step 705, then the value of that pixel the maximnm value. Thus, a determination is made of 

b'C j) is set equal to the max value. After step 710 or if whether b'(i j) is greater than the iha x i m nm value, 

a negative determination is made at step 706, then a If an affirmative determination is made at step 906, 

detenninatioa is made at step 70S of whether the value then for each pixel in the array b^ij) which is greater 

of the pixel of the modified blue noise mask b'G J) is less 5 than the maximum value, its value is set to the ni a ximum 

than the minimmn value. If an affirmative determination value at step 910. 

is made at step 708, then the value of that pixel b'(i j) is I^ a negative determination is made at step 906 or 

set equal to the minimum value at step 712. After step after step 910, a determination is made at step 908 of 

712 or if a negative determination is made at step 708, which pixel value, if any, in the modified blue noise 

the modified blue noise mask is assembled by writing 10 mask b^ij) is less than the niinimum value. If an affirma- 

the values of such modified blue noise mask b'(ij) into ^ive determination is made at step 908, then for each 

memory at step 714. pi^^^ value in the modified blue noise mask b'Cij) which 

FIG. 8 shows a diagram of the flow chart for generat- ^ less than the minimum value, that value is set equal to 

ing a modified high contrast version of the blue noise minimum value at step 912. ... 

mask. At step 800, the values of the blue noise mask 15 After step 912 or if a negative determmation is made 

b(i,j)=g are read. At step 802, the constants a, b, c values of the modified blue noise mask 

which are to be utilized in step 804 are obtained from ^'Oj) ^® memory at step 914. 

the user, as weU as the max and min values. The con- instant method can also be apphed to color half. 

stants a, b and c are chosen so as to generate unique ^^6' mdependently thr^holdmg each one of the 

^mathematical mapping. The center point of inflection is ^ <^omponent colors agamst the blue noise mask and then 

given by the coiitant a, usually chosen near 128 for a overprmtmg. In tl^t manner, die blue noise mask dis- 

« Kit «,oci. Ti,^ c*lo«««o «f +1,0 z.™,** « closed herem can be apphed simply to die component 

l^t^tw^..^^^ Et of Se ™ ff <^lors of ROB, CMYK and others for color printing, 

given by the constant c and the offset of the curve, if ^^^^^ ^ ^^^^ ^^^^ 

d^ed. by the constant b- 25 scanner 400 shown in FIG. 4 to input three separate files 

Direcdy mappmg occjus at step 804 m which each g.^.^ P ^^^^ ^een com- 

pixel value for the modified blue no^se mask ^ (ij) is ^^^^ p^^^ displaying such image, a 

calculated to be equal to (bCij)-a)yc3+b. This func ^^^^ generated accordkg to the method 

Hon replaces a hnear mputHDutput relation with a steep- disclosed herein can be applied separately to each of the 

cned. non-hnear relation. The hnear mput-output rela- 30 ^^^^ ^ ^^^^ ^ 

tion being replaced is shown, for example, in FIG. 10 displayed on a binary RGB video screen or 

which shows the CDF versus the value of pixels for a printed 

linear blue noise mask. FIG. 11 shows the CDF versus discovered that an improvement in 

the value of pixels for a non-linear, high contrast blue ^y^^y ^f^^ displayed RGB image can be achieved 

noise mask, after applying the CDSC direct mapping 35 when the pixels ofthe blue noise mask are shifted by one 

with punch. pixel when used on the different color planes. For ex- 

A determination is then made at step 806 of whether ample, the (ij) pixel of the blue noise mask can be 

each new pixel value b'(ij) is greater than the maximum shifted such that the blue noise mask which is applied to 

value. If an affirmative determination is made at step i^g^ color plane has each of its pixel values 

806, then the new pixel value b'(ij) is set equal to the 40 shifted up by one pixel when halftoning the red image, 

maximum value at step 810. After step 810 or if a npga- halftoning the blue image, the blue noise mask 

tive determination is made at step 806, then a determina- h^s each of its pixel values shifted by one pixel to the 

tion is made at step 808 of whether the new pixel value side. In that manner, the color energy is spread out over 

b'(i,j) is less than the miniinum value. a larger space. It should be noted, however, that no 

If an affirmative determination is made at step 808, 45 shifting of the RGB images with respect to each other 

then the new pixel value b'(ij) is set equal to the mini- occurs, as that would cause blurring of the resulting 

mum value. After step 812 or if a negative determina- image. Instead, as discussed above, the respective blue 

tion is made at step 808, then the values ofthe modified noise masks used to halftone the red and blue images 

blue noise mask b'O j) are written mto memory at step have their pixel values shifted by one pixel up or to the 

814. 50 side, respectively. 

FIG. 9 shows a diagram of the flow chart for the xhe foregoing method is shown in FIG. 12 which is 

modification of a blue noise mask to produce a CUSC a dxa^am of a flow chart for implementing halftoning 

low contrast version of the original blue noise mask. of a color image using the foregoing method The color 

The values of the blue noise mask b(ij)=g are read at image 1200 to be scanned is scanned by a scanner 400 to 

step 900 and then the constants a, b and c are obtained, 55 produce three contniuous tone color planes 1202, 1204 

as well as the maximum and miniTnum values at step and 1206, one for each of the three colors, red, green 

902. Such constants and maximum and minimum values and blue, respectively. The values in each of these 

are provided by the user. planes or arrays are added to the blue noise mask as 

At step 904, a direct mapping process is accomplished shown in FIG. 6. Alternatively, a comparator can be 

in which the array of values of the modified blue noise 60 used, as shown in FIG. 4. However, the (i J) pixel of the 

mask b'(ij) is calculated as equal to cbrt (b(i,j)— b)*- red image is compared against the shifted (i-f-l j) pixel 

c-ha. This ftmction, where "cbrt" stands for the cube ofthe blue noise mask. The CkJ) pi^^l ofthe green image 

root, changes a linear input-output relation to a low is compared against the (iJ) pixel ofthe blue noise mask, 

contrast, non-linear input-output relation. The constants The (i j) pixel of the blue image is compared against the 

a, b and c give the of&et, center point and gain, respec- 65 OJ + 0 pi^^ c>f the blue noise mask. Such comparisons 

tively. take place at steps 1208, 1210 and 1212. 

At step 906, a determination is made for each pixel of Each of these three planes or arrays are thresholded 

the blue noise mask of whether its value is greater than or are printed on high gamma film at step 504 in order 
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to produce three halftone images at steps 1214^ 1216 and 
1218. The halftone image h/jCu) represents the (i j) pixel 
of the halftone red image. The elements at 1216 and 
1218 likewise represent the halftone of the green and 
blue images, respectively. These images are combined 5 
at step 404 by a printer or display in order to produce a 
three color halftone image at step 1220. It should be 
noted that the halftoning of each of the components of 
the color image uses a shifted version of the blue noise 
mask as disclosed herein to produce a useful spreading 10 
out of the color without introducing blur into the half- 
tone image. 

A variation of this method of separately halftoning 
the different color planes of an image is to use the in- 
verse of a blue noise mask for one color, where the IS 
inverse is defined as (Maxval-BNM) for each pixel The 
Maxval is the maximum value of the blue noise mask 
(for example, 255 for the 8-bit blue noise mask). To 
produce the modified blue noise mask, the value of each 
pixel is subtracted &om the Maxval to produce each 20 
new value for each pixel. This inverse process can be 
thought of as exchanging the "peaks and valleys of the . 
blue nobe mask," and results in the spreading out of the 
energy in the color pattern of the color image. 

Many modifications of this process are likewise possi- 25 
ble, such as different one shift patterns, or multiple shift 
patterns of the pixel values of a generated blue noise 
mask. While it has been found, for example, in low 
resolution 300 dpi systems, that one pixel shHt produces 
a visually pleasing result, ten shifts have been found to 30 
produce an annoying correlation between different 
color dots while larger shifts of about forty-five pixels 
produce an acceptable, uncorrelated pattern of colored 
dots. 

Of course, other possible modifications can also pro- 35 
duce acceptable results, for example, the placing of one 
color plane blue noise mask at 45" with respect to an- 
other as is done in conventional four color printing. 
However, it has been found that the simple one-shift 
pattern described above is most effective at producing a 40 
pleasing, isotropic, non-clumpy moire resistant pattern 
with some spreading out of the color or ink but with no 
blurring of the image. Similar principles apply of course 
for CMYK color printing, where the halftoning can be 
done on the C, M and Y color images, and then the 45 
black image (K) can be added as necessary according to 
conventional color printing models. 

Although only a preferred embodiment is specifically 
illustrated and described herein, it will be readily appre- 
ciated that many modifications and variations of the SO 
present invention are possible in light of the above 
teachings and within the purview of the appended 
claims without departing from the spirit and intended 
scope of the invention. 

What is claimed is: S5 

1. A method for the halftoning of color images, com- 
prising the steps of utilizing, in turn, a pixel-by-pixel 
comparison of each of a plurality of color planes of said 
color image against a blue noise mask in which the blue 
noise mask is comprised of a random non-deterministic, 60 
non-white noise single valued function which is de- 
signed to provide visually pleasing dot profiles when 
thresholded at any level of said color images, wherein a 
plurality of blue noise masks are separately utilized to 
perform said pixel-by-pixel comparison and in which at 65 
least one of said blue noise masks has its pixels shifted by 
at least one pixel prior to performing said pixel-by-pixel 
con^)arison. 



2. The method of claim 1, wherein said blue noise 
mask has a maximum size value and wherein a plurality 
of blue noise masks are separately utilized to perform 
said pixel-by-pixel comparison and in which at least one 
of said blue noise masks has its pixels modified by re- 
placing the value of each pixel with a value correspond- 
ing to the difference between said fnaximnm size value 
and the pixel value prior to performing said pixel-by- 
pixel comparison. 

3. A method of halftoning a color image by utilizmg 
a pixd-by-pixel comparison of said color image against 
a blue noise mask array stored in a first memory, com- 
prising the steps of: 

a) scanning said color image to be halftoned to create 
a plurality of color plane image arrays, said scan- 
ning being performed on a pixd-by-pixel basis; 

b) storing each of said plurality of color plane image 
arrays in a second memory; 

c) comparing, on a pixel-by-pixel basis, the value of 
each corresponding pixel m said blue noise mask 
array and each of said plurality of color plane 
image arrays to produce a plurality of binary color 
plane image arrays; 

d) converting said plurality of binary color plane 
image arrays to the desired halftoned color image; 
and 

e) producing at least one modified blue noise mask 
array in addition to said blue noise mask array prior 
to performing said comparing step; said at least one 
modified blue noise mask array being generated by 
shifdng the pixels of said blue noise mask array by 
at least one pixel position from their position in said 
blue noise mask array. 

4. The method of claim 3, wherein said at least one 
modified blue noise mask array is compared, on a pixel- 
by-pixel basis, to the value of each corresponding pixel 
in at least one of said plurality of color plane image 
arrays. 

5. The method of claim 3, further including the step 
of producing at least one modified blue noise mask array 
in addition to said blue noise mask array prior to per- 
forming said comparing step; said at least one modified 
blue noise mask array having of a predetermined maxi- 
mum value and being generated by subtracting the 
value of each pixel of said blue noise mask array from 
said maximum value to produce said modified blue 
noise mask array. 

6. The method of claim 3, wherein said at least one 
modified blue noise mask array is independent and un- 
correlated with any other blue noise mask array. 

7. A method for generating a modified blue noise 
mask array in which said modified blue noise mask 
array may be used to generate a more pleasing halftoned 
image than in a non-modified blue noise mask array, 
comprismg the steps of: 

a) reading each of the values of the blue noise mask 
array to be modified; 

b) receiving maximum and minimimi predetermined 
values input by a user; 

c) modifying each value of said blue noise mask array 
using at least one of said maximum and miniTniiTn 
predetermined values and a direct value mapping 
function; 

d) comparing each of said values of said modified 
blue noise mask array to said maximum predeter- 
mined value and replacing said value if it exceeds 
said maximum predetermined value; and 
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e) comparing each of said values of said modified blue 
noise mask array to said minimum predetermined 
value and replacing said value if it is less than said 
minimum predetermined value. ^ 

8. The method of claim 7, further including the step 
of storing said modified blue noise mask in a memory* 

9. The method of claim further including the step 
of using said stored modified blue noise mask to halftone 
an image. 

10. The method of claim 7, wherein said maximum 
and minimum predetermined values input by said user 
and said direct value mapping fimction are selected to 
compensate for characteristics of at least one of printer 15 



10 



20 



and display devices used by said user so as to form a 
more pleasing halftoned image. 

11. A method for the halftoning of color images, 
comprising the steps of utilizing, in turn, a pixel^y-pixel 
comparison of each of a plurality of color planes of said 
color image against a blue noise mask in which the blue 
noise mask is comprised of a random non-deterministic, 
non-white noise single valued function which is de- 
signed to provide visually pleasing dot profiles when 
thresholded at any level of said color images, wherein a 
plurality of blue noise masks are separately utilized to 
perform said pixel-by-pixel comparison and in which at 
least one of said blue noise masks is independent and 
uncorrelated with the other blue noise masks. 
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